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^STRACT 


Analytical  solution  of  the  approximation  problem  of  network  syn- 
thesis entails  a large  amount  of  computational  work.  This  report  describes 
an  experimental  method  of  solution  by  means  of  a computer  based  on  the 
electrostatic  analogy  in  which  voltages  in  a conducting  sheet,  set  up  by 
line  currents,  represent  network  immittance  functions.  The  resulting  approxi- 
mating function  is  described  in  terms  of  its  pole  and  zero  locations  and  a 
multiplicative  constant.  The  method  of  operation  of  the  analog  computer  is 
simple.  The  poles  and  zeros  are  moved  around  in  the  conducting  medium  until 
the  desired  immittance  is  obtained.  The  continuous  scan  feature  of  the  com- 
puter  makes  possible  visualization  of  the  response  and,  what  is  more  important, 
immediately  shows  the  effect  of  moving  any  of  the  poles  or  zeros.  A conver- 
gence procedure  for  the  perturbation  of  these  critical  points  has  been  de- 
veloped so  that  an  apprcxi mat-ion  within  the  required- accuracy  is  often  possible 
within  a very  short  period  of  time. 

Approximation  from  the  viewpoint  of  economization  is  carried  out 
with  the  objective  of  reducing  the  number  of  critical  points  of  a known 
approximating  function  in  order  to  reduce  the  number  of  required-  elements  in 
the  network.  The  process  entails  the  removal  of  a pole  or  zero  and  redistri- 
bution of  the  remaining  critical  points  to  return,  as  closely  as  possible,  to 
the  required  response.  The  removal  of  one,  or  a pair,  of  critical  points, 
perturbation  of  the  remaining  ones,  and  adjustment  of  the  multiplicative  con- 
stant is  then  repeated.  The  process  is  stopped  just  before  the  number  of  poles 
and  zeros  has  been  so  greatly  reduced  as  to  make  the  approximation  within 
allowed  deviation  unobtainable. 

Positive  real  functions  are  guaranteed  to  give  physically  realizable 
driving-point-  immittances.  To  insure  PRF  results,  the  analog  approximation  is 
made  to  the -real  part  of  the  driving-point  immittance.  The  computer  is  well- 
suited  for  approximations  to  the  real  or  imaginary  parts,  as  well  as  the  ampli- 
tude function. 

Approximations  were  carried  out  to  the  transcendental  function  aris- 
ing from  the  input  admittance  of  a shorted  transmission  line.  The  synthesized 
networks  were  found  to  give  similar  deviations  than  known  partial  fractions 
expansions.  Furthermore,  the  analog  gave  networks  with  fewer  elements. 
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J c In  t r o due tic  n 

Although  t.he  'lieory  of  synthesis  of  passive,  linear,  lumped-parameter 
networks  is  well-known,  there  is,  as  yet,,  no  completely  satisfactory  solution 
to  the  approximation  problem.  Usually,  a synthesis  can  be  carried  out  after 
the  immittance*  function  can  be  described  by  means  of  an  algebraic  expression 
of  a specified  form.  Often,  the  graphical  representation  of  the  function  is 
given  and  the  problem  then  becomes  one  of  finding  the  equation  of  the  curve,  in 
the  required  form,  which  will  approximate  the  desired  function  to  a required 
degree  of  accuracy.  Analytical  methods  of  solution  exist,  but  the  amount  of 
computational  work  is  generally  very  large. 

This  report  describes  an  analog  computer  which  provides  an  experimen- 
tal means  of  solving  the  approximation  problem.  The  procedure  is  essentially 
ore  of  trial- and-err or,  but  converges  very  rapidly.  The  solution  is  found 
directly  in  terms  of  the  locations  of  the  poles  and  zeros  of  the  function  re- 
presenting the  required  response. 

The  methods  used  here  to  solve  the  approximation  problem  are  made 
possible  because  of  an  analogy  which  exists  between  electrostatics  and  the 
type  of  rational  analytical  functions  which  arise  in  ohe  network  theory  (Refer- 
ences 1,  Uj  and  8).  Thus,  the  highly  developed  mechanism  of  potential  theory 
becomes  available  for  the  solution  of  circuit  problems.  The  common  and  easily 
developed  "feeling*  for  electrostatic  problems  may  be  applied  to  give  useful 
qualitative  results  in  network  synthesis. 


II.  The  Electrostatic  Analogy 

Any  immittance  function  can  be  described  in  terms  of  its  poles  and 
zeros  and  a multiplicative  constant, 


F(P) 


H(p  - P2)  (P  - Pj^)  ••• 
(P  ~ Px)  (P  - P3)  ••• 


(1) 


where, 

p - '-r  + j<&  is  the  complex  frequency  variable 

H is  the  multiplicative  constant 

P2,  » P?n  are  the  zerios  F(p) 

p-j,  py  ...  . P2n-1  are  the  Poles  of  f(p) 

The  logarithm  of  the  absolute  value  of  the  immittance  is 

* In  this  report,  the  general  term  "immittance  function"  includes  impedance  or 
admittance,  of  both  driving-point  or  transfer  type. 
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ln|Ftp)|  - In  H ln|p  p?n| ‘ ^InlP'Pjn I (2) 

and  is  called  the  gain  function  of  the  network.  Note  that  the  gain  function, 
at  any  point  p..  is  equal  to  the  sum  of  the  logarithm  of  the  distances  from  p 
to  the  zeros,  minus  the  sum  of  the  logarithm  of  the  distances  to  the  poles, 
plus  a constant. 

Now,  let  us  examine  the  voltages  in  a uniform,  infinite,  conducting 
sheet  pierced  by  line  current  sources. 

Consider,  a finite  region  of  the  infinite  conducting  plane  of  thick- 
ness o , as  shown  in  Fig.  MRI-13292»  0 is  the  origin,  where  the  current  I is 
introduced  into  the  plane.  Q is  any  point  whose  distance  from  the  origin  is  q, 
and  K is  a reference  point  located  at  distance  k from  0„ 

Since  the  conducting  sheet  is  infinite,  the  current  emanates  radially 
from  the  source  at  0,  But, 


JA 


(3) 


The  assumption  is  made  that  there  is  no  current  flow  in  the  S direction.  Also, 


7/ 


(h) 


hence . 


IR 


J j>  1 


f 


(5) 


where,  J * current  density 

^ * resistivity  of  paper 
1 - length  of  conducting  path 
A * area  presented  to  the  current 

Then,  at  any  radial  distance  r,  the  change  in  voltage  dV  for  a 
differential  change  in  radius  dr  i s 


dV 


1 ^ dr 

2 it  i r 


(6) 


Integrating  from  K to  Q results  in 

V *.  In  »> 

KQ  2n£  m q 

” C - .^2,  In  q 

?.TI 


(7) 
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where  Ro  is  the  resistance  per  square  of  the  paper  and  C is  a constant. 

It  is  seen,  then,  that  for  a constant  current  source,  the  voltage 
at  any  point  is  proportional  to  the  logarithm  of  the  distance  to  the  current 
source. 


If  now,  the  conducting  sheet  represents  the  p-plane,  and  if  posi- 
tive currents  are  introduced  at  points  corresponding  to  the  poles  and  negative 
currents  are  introduced  at  points  corresponding  to  the  zeros,  all  currents 
being  of  the  same  magnitude,  the  resultant  voltage  becomes,  by  superposition, 

(8) 


IR 


A + 


o 

Tn 


- ■|Vlp-p2n_1l 


Eq.  (8)  is  Identical  in  form  with  Eq.  (2)  for  the  gain  function  of 
a network.  The  constant  multiplier,  IR_/2jt,  determines  the  scale  factor. 

The  constant  A determines  the  zero  level  of  the  voltage,  i.e.,  the  impedance 
level  of  the  gain  function. 


III.  The  Logarithmic  Transformation 

Double-layer,  circular,  electrolytic  tanks  representing  the  p-plane 
have  been  used  (References  2 and  3)  as  a tool  in  network  analysis.  The  use 
of  carbon  impregnated  paper  as  a conducting  medium  is  facilitated  by  the 
logarithmic  transformation  (Reference  6) 


W *=  In  p 

- ln|p{  + ,1$ 

- U + j V 


where,  p 


be 


.1* 


cr+  > 


(9) 


The  Riemann  surfaces  of  the  p~plane  are  all  represented  by  the  in- 
finite W-plane  sheet.  The  single  Riemann  surface 


01':  $ <2n 


made  by  taking  the  branch  cut  in  the  p-plane  along  the  positive  real  axis, 
corresponds  to  a semi- in finite  Vf-plane  strip  of  width  2n.  The  transformation 
is  shown  in  Fig.  MRI-13291*. 
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The  advantages  of  the  logarithmic  transformation  led  to  the  con- 
struction of  the  physical  analog  in  the  logarithmic  plane.  These  advantages 
are: 

1.  The  logarithmic  coordinates  are  convenient  for  many  problems. 

2.  A sizeable  range  of  frequencies  can  be  represented  in  a convenient 
size  sheet. 

3-  Calibration  is  simple  using  known,  constant- slope  response  curves. 

U.  The  accuracy  is  uniform  over  the  useful  region  of  the  plane. 

5.  Errors  due  to  the  finite  size  of  the  plane  are  easily  made  negligi- 

ble . 

The  main  advantages  of  using  the  conducting  paper,  rather  than  an 
electrolyte.  Are? 

1.  The  paper  ia  net  affected  by  shock,  whereas  disturbances  of  the 
electrolyte  cause  ripples  and  consequent  distortion  of  the  field. 

2.  Use  of  direct  current  is  possible  since  there  is  no  dissociation 
problem. 

3-  Paper  is  more  conveniently  handled  then  a fluid 


IV.  Symmetry  Conditions 

A.  Symmetry  abcut  the  Real  Axis 

Since  network  iimittances  are  analytic  functions  expressable  as 
ratios  of  polynomials  in  p with  real  coefficients,  poles  and  zeros  appear 
either  as  complex  conjugate  pairs  or  on  the  real  axis.  This  means  that  the 
W- plane  will  have  current  sources  and  sinks  placed  symmetrically  about  the 
n(real)  axis.  Thus,  nc  current  will  flow  across  the  Tt-axis  and  the  plane  may 
be  cut  along  this  line  without  disturbing  the  field,  or  voltage  distribution. 
Only  one  half  of  the  transformed  plane  is  then  used  and  it  contains  all  the 
necessary  information.  Poles  And  zeros  on  the  real  axis  are  reduced  to  half 
their  normal  value  since  only  half  the  current  would  normally  flow  into  the 
half- plane  used  in  the  analog. 

This  symmetry  condition  is  shown  in  Fig.  MRI-13295-A. 
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B.  Symmetry  About  the  Imaginary  Axis 


Normally,  the  poles  and  zeros  of  an  immittance  are  not  symmetrically 
placed  about  the  imaginary  axis.  If  the  symmetrical  critical  points  are  in- 
troduced into  the  analog,  the  voltage  produced  along  the  n/2  (imaginary)  axis, 
that  is,  the  gain,  is  simply  doubled  (Reference  $).  Thi3  result  can  be  proven 
mathematically  in  the  following  fashion  (Reference  U) : Let  the  original  im- 

mittance (normalized)  be 


F(p) 


*(P-P2n) 


(10) 


as  in  Eq.  (1).  The  logarithm  of  the  magnitude  of  the  immittance  is  then  the 
voltage  measured  along  the  imaginary  axis  or. 


>~p2n 

•>*p2»-l 


With  the  addition  of  the  image 

Vp)  ‘ 


critical  points,  F(p)  is  changed  to 
*(P-P2c)  (P-P^) 

Jl(p’p2m-1)  (P-Pan-l* 


(11) 


(12) 


Now  the  logarithm  of  the  magnitude  of  this  function  along  the  imaginary  axis  is 
the  measured  voltage 


*-p2n 

j“-p2m-l 


>-p2n 

'l!fcp2m.l 


(13) 


The  remit  shows  that  the  magnitude  measurement  is  simply  doubled. 

The  imposed  symmetry  about  the  imaginary  axis  makes  it  possible  to 
cut  the  analog  plane  in  half  again,  as  did  the  symmetry  about  the  real  axis. 
This  time  critical  points  on  the  imaginary  axis  are  reduced  to  half-currents 
and  only  one  quarter  of  the  original  p-plane  is  represented  in  the  analog. 
This  symmetry  condition  is  shown  in  Fig.  MRI=1329.5“B.  Furthermore,  the 
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imaginary  axis  now  represents  a stream  line  (i.e.,  no  current  flows  across 
it).  Since  all  current  flow  lines  approach  the  imaginary  axis  assymptotically, 
and  since  current  and  voltage  are  conjugate  functions  in  the  analog,  then  the 
constant  voltage  lines  are  all  perpendicular  to  the  imaginary  axis,  which 
means  that  the  voltage  measuring  probe  does  not  have  to  be  so  accurately  posi- 
tioned. 


V.  Description  of  the  Analog  Computer 

The  computer  is  mounted  on  a standard  size  machinist's  table  having 
three  shelves  which,  with  the  addition  of  an  oscilloscope,  makes  a complete, 
easily  portable  unit.  A drawing  board  mounted  on  the  top  shelf  acts  as  the 
support  for  the  conducting-paper  plane.  The  remaining  components  are  placed 
on  the  other  two  shelves. 

Currents  are  introduced  into  the  conducting  plane  by  means  of  cur- 
rent probes.  These  probes  are  3imply  steel  needles  pressed  into  brass  sleeves. 
The  sleeves  facilitate  the  attachment  of  leads.  The  current  probes  are  fed 
from  d-c  sources  which  are  constant  voltage  supplies  of  either  plus  or  minus 
150  volts.  1$0  kilohms  is  inserted  in  series  with  each  probe  for  full  currents 
and  300  kilohms  for  half- currents.  Thus,  the  full  currents  are  of  the  order  of 
one  ndlliampere.  This  was  chosen  as  a convenient  value  since  it  leads  to  con- 
venient values  of  scanning  voltage.  In  calibrating  (see  p.  10),  the  voltage 
gradient  along  the  real  frequency  axis  is  approximately  0.6  volts  per  decade 
of  frequency  and  is  readily  handled  by  the  display  system.  Lower  current 
sources  would  require  further  voltage  amplification  while  higher  currents  tend 
to  heat,  and  possibly  burn,  the  conducting  paper.  Series  potentiometers  of 
20  kilohms  are  used  to  adjust  all  currents  to  the  required  values.  These  are 
necessary  since  the  impedance  seen  by  a current  probe  depends  on  its  position 
in  the  conducting  sheet  as  well  as  on  the  total  number  of  probes  in  the  sheet. 
In  order  to  preserve  the  drawing  board,  a replaceable  cork  sheet  is  placed 
beneath  the  conducting  paper.  The  current  probes  are  held  in  place  by  pressing 
into  the  cork  sheet. 

The  voltage  along  the  real  frequency  axis  is  scanned  by  a rolling 
voltage  probe.  The  voltage  probe  comprises  a ball  bearing  with  a copper  sleeve 
pressed  on  to  the  outside  of  the  race.  The  sleeve  is  circular  in  cross-section 
and  hence,  point- con tact  is  made  with  the  paper.  The  voltage  probe  carriage 
is  supported  on  a length  of  rectangular  waveguide  and  is  driven  across  the 
sheet  by  means  of  a motor  and  belt  drive.  Reversing  micro- switches  are  mounted 
at  each  end  of  the  guide.  The  driving  motor  is  a commutator  a-c  type  with  a 
forced  air  cooling  system.  Since  the  total  scanning  time  of  the  paper  is  short, 
a cyclic  duty  motor,  instantly  reversible  every  five  seconds,  is  required.  A 
high  torque  to  inertia  ratio  and  a cooling  fan  prevent  overheating. 

The  voltage  at  the  pick-up  (scanning)  probe  is  displayed  on  a Dumont 
30l*H  cathode  ray  oscilloscope.  The  scope's  d°c  amplifiers  have  sufficient  gain 
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for  the  voltages  encountered  so  that  no  pre-amplifiers  are  necessary.  Horizon- 
tal deflection  for  the  display  system  is  obtained  from  a Helipot  precision 
linear  potentiometer  driven  from  the  voltage  probe  carriage  by  a continuous 
belt  drive.  A tension  adjusting  coil  spring  is  inserted  in  the  belt  to  mini- 
mize backlash.  Drag  on  the  probe  carriage  ha3  been  reduced  to  a minimum  es- 
tablished by  the  friction  of  the  pulleys  and  potentiometer  bearings.  Since  the 
output  from  the  analog  plots  a logarithmic  gain  function  gainst  the  logarithm 
of  frequency,  the  face  of  the  cathode  ray  tube  has  been  calibrated  in  log- log 
coordinates.  This  was  done  by  spraying  the  tube  face  with  an  acrylic  coating. 

A log-log  grid  was  then  inscribed  with  a needle-point  on  the  surface  thus 
formed.  The  scratch  marks  were  filled  with  marking  crayon  to  make  permanent 
coordinates  directly  on  the  face  of  the  tube.  A long-persistance  tube  is  used 
in  order  to  retain  the  trace  for  as  long  a time  as  possible.  This  makes  the 
effect  of  movement  of  critical  points  readily  observable.  The  computer  is  thus 
a continuous  scan  mechanism.  A picture  of  the  analog  computer  (without  the 
display  oscilloscope)  is  shown  in  Fig.  MRI-13296. 


VI-  Errors  in  the  Cc^npruter 

The  percentage  error  depends,  somewhat,  on  the  particular  problom 
involved.  There  are  certain  errors,  however,  which  are  inherent  to  the  device. 
These  are  analyzed  below.  The  principle  sources  of  error  are: 

1.  Nonuniformity  of  the  conducting  paper 

2 „ Paper  shift 

3 Finite  size  of  the  plane 

U.  The  probes 

5.  The  oscilloscope 

a.  nonlinearity  of  the  amplifiers 

b.  parallax 

c . pickup 

1 . __  Nonuniformity  of  the  Conducting  Paper 

Changes  in  the  resistivity  of  the  paper  may  be  caused  by  a change 
in  the  spacing  of  the  rollers  while  the  paper  is  being  made  or  by  impurities 
and  inhoraogeneities  in  the  material.  Errors  due  to  a change  in  roller  spacing 
are  completely  negligible  since  the  long  dimension  of  the  conducting  sheet 
used  in  the  analog  is  taken  from  the  width  of  the  roll  and,  usually,  only  3»Ul 
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inches  of  length  along  the  roll  are  needed.  The  latter  cause  of  nonuniformity 
was  investigated  by  measuring  the  resistance  per  square  of  the  conducting 
sheet.,  A pair  of  coaxial  circular  conductors  were  placed  in  contact  with  the 
conducting  medium.  The  resistance  between  these  conductors  is  directly  pro- 
portional to  the  resistance  per  square  of  the  conducting  sheet,  the  propor- 
tionality constant  being  determined  by  the  geometry.  The  area  between  the  con- 
ductors is  approximately  1.3  square  inches.  Hence,  this  method  is  a measure  of 
the  average  resistivity  over  a comparatively  small  area.  Tests  were  carried 
out  on  a 5 foot  sheet  cut  from  the  roll,  which  is  about  31  inches  wide.  About 
200  readings  were  taken,  or  one  reading  every  three  inches.  The  reailts  showed 
that  the  paper  has  a resistance  of  1800  ohms  per  square  with  variations  of  less 
than  0.U  percent. 


2 . Paper  Shi  ft 

Paper  shift  was  an  extremely  troublesome  problem.  In  order  to  assure 
good  contact  of  the  rolling  voltage  probe,  it  is  necessary  to  spring  load  the 
contacting  wheel.  Unless  the  conducting  paper  is  securely  held,  the  voltage 
probe  will  shift  the  paper  as  it  moves  across  the  plane.  The  spot  on  the 
cathode  ray  tube  then  does  not  retrace  the  same  curve  when  the  voltage  probe 
scans  in  the  opposite  direction.  Paper  shift  has  been  eliminated  with  the  use 
of  spring-fingers,  at  both  ends  of  the  conducting  sheet,  arranged  sc  as  to 
flatten  the  sheet  and  put  it  into  tension. 

3°  Finite  Size  of  the  Plane  - 

In  the  logarithmic  plana  no  errors  are  introduced  by  the  finite 
width  of  the  paper,  providing  that  the  ratio  of  width  to  length  is  such  that 
the  conformality  of  the  transformation  is  preserved.  The  conducting  sheet 
was  chosen  20  inches  long  as  a convenient  size.  Four  decades  of  frequency  are 
used,  making  each  decade  five  inches  long.  The  corresponding  width  must  re- 
present a change  of  2n  radians.  Hence,  for  the  width  x. 


- J _ 

2*  lnlO 


(Hi) 


Therefore,  x is  13.65  inches.  But  the  analog  represents  only  a quarter  of  the 
p-plane.  In  practice,  then,  the  conducting  sheet  is  3»l4l  inches  wide. 

The  transformation  requires  the  use  of  a send -infinite  sheet.  There 
is  an  error  introduced  in  using  a finite  length  to  represent  the  logarithmic 
plane.  Consider  the  field  far  removed  from  a finite  region  where  all  the  poles 
and  zeros  are  located.  The  return  for  any  excess  current  introduced  into  the 
plane  is  taken  at  infinity.  Then  the  equipotentials  far  from  the  furthest  es- 
sential point  are  constant  amplitude  lines  or  lines  perpendicular  to  the  length 
of  the  conducting  sheet.  Equipotentials  are  painted  across  the  sheet  with 
silver  paint  to  represent,  at  the  two  ends  of  the  paper  respectively,  the  origin 
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and  the  point  at  infinity.  Analysis  shows  (Reference  h)  that  if  the  equi- 
potential  representing  the  point  at  infinity  is  placed  at  seven  times  the 
distance  to  the  furthest  singular  point,  the  errors  due  to  the  finite  size 
of  the  plane  will  never  be  greater  than  one  percent.  In  the  computer,  this 
ratio  has  been  increased  to  ten  as  an  added  precaution.  The  same  situation 
occurs  at  the  other  end  of  the  strip,  and  it  is  likewise  satisfied  by  leaving 
an  unused  decade  at  the  zero  end.  The  sheet  is  four  decades  long.  Thus,  the 
two  center  decades  are  available  for  pole  and  zero  location,  or  the  frequency 
range  is  100  to  one. 

k*  The  Probe s 

Because  the  current  probes  are  finite  in  size  they  distort  the  field 
by  disrupting  the  homogeneity  of  the  conducting  medium.  Holes  left  in  the 
paper  after  the  removal,  of  the  neodle  probes  have  the  same  effect.  The  die- 
tortion  of  the  field,  however,  is  negligible  at  distance  greater  than  ten  times 
the  radius  of  the  probe  (Reference  U).  Since  the  probe  radiuB  is  .0125  inches, 
distortion  due  to  prcbes  is  not  noticeable  at  distances  of  more  than  one -eighth 
inch.  A similar  condition  holds  for  the  voltage  probe  since  the  rolling  wheel 
does  not  make  actual  point  contact  but  has  approximately  a circular  contact  area 
of  small  radius.  A further  error  is  Introduced  if  the  voltage  probe  does  not 
accurately  follow  the  line  of  the  real  frequency  axis.  This  error  has  been 
minimized  in  the  computer  by  making  the  constant  voltage  lines  perpendicular 
to  the  scanning  axis  (p.  6). 

5.  The  Oscilloscope 

Nonlinearity  of  the  amplifiers  can  be  kept  to  a minimum  of  3 percent 
by  proper  adjustment  of  the  scope  and  use  of  fresh  tubes.  The  center  three 
inches  of  the  five  inch  tube  used  is  almost  perfectly  linear;  hence,'  non- 
linearity is  not  an  important  error  in  the  center  of  the  tube. 

Parallax  error  was  minimized  by  inscribing  the  log-log  coordinates 
directly  on  the  face  of  the  cathode  ray  tube  (p.  7).  The  separation  of  the 
spot  and  grid  is  then  only  the  thickness  of  the  glass  face  of  the  tube. 

Pickup  produces  error  since  it  tends  to  give  thickness  to  the  spot, 
making  fine  focusing,  and  hence  exact  location  of  the  point,  impossible. 

In  roost  cases,  the  largest  errors  are  introduced  by  the  display 
system,  i.e.,  the  oscilloscope.  Because  of  diffusion,  nonlinearity  and  parallax 
the  spot  cannot  be  located  closer  than  within  one-sixteenth  of  an  inch.  This 
introduces  an  error  on  the  log-log  grid  of  no  more  than  five  percent.  The 
maximum  overall  error  is  then  approximately  six  percent. 

If  the  function  in  question  has  a critical  point  (zero  or  pole) 
within  an  eighth  of  an  inch  of  the  real  frequency  axis,  the  display  is  not 
valid  immediately  within  the  vicinity  of  the  current  probe.  This  can  be 
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readily  remedied,  however,  by  calculating  a point  or  two  on  the  curve  in 
question. 


VII.  Method  of  Operation 

The  method  of  operation  of  the  computer  is  comparatively  simple. 

The  desired  form  of  the  gain  function  is  drawn  directly  on  the  face  of  the- 
scope  with  China  marking  crayon.  The  poles  and  zeros  are  then  moved  around 
in  the  conducting  sheet  until  the  picture  on  the  cathode  ray  tube  agrees  with 
the  desired  one  within  specified  tolerance. 

In  this  manner,  the  approximation  problem  of  network  synthesis  can 
be  solved  in  comparatively  a very  short  time.  In  many  cases,  the  correspond- 
ing analytical  procedure  would  require  several  days.  Furthermore,  problems 
requiring  a large  number  of  essential  points  are  handled  with  almost  the  same 
facility  as  the  simple  problems.  This  is  not  at  all  true  of  the  analytical 
approach. 


VIII.  Calibration  Procedure 

Calibration  of  the  log-amplitude  scale  presents  no  difficulty.  Es- 
sentially the  procedure  involves  setting  up  a known  amplitude  function  in  the 
analog  and  setting  the  gain  controls  to  see  the  proper  response  on  the  cathode 
ray  tube.  The  simplest  response  is  that  due  to  one  pole  (or  one  zero).  If 
the  pole  is  located  at  the  origin,  the  response  is  a constant  slope  line. 

That  is, 

F(«)  - - (15) 

00 

From  which, 

slope  - d(log_[FjJ-  . _ ! (16) 

d(log  o>) 

The  procedure,  then,  is  to  set  the  horizontal  positioning  and  gain  controls  so 
that  the  position  of  the  voltage  probe  corresponds  to  the  frequency  scale, 
then  scanning  the  output  voltage  from  a simple  pole  at  the  origin  and  setting  the 
vertical  gain  control  so  as  to  see  a unit  slope  line.  For  convenience,  this 
slope  has  been  permanently  inscribed  on  the  face  of  the  cathode  ray  tube.  It 
should  be  noted  that  a simple  pole  at  the  origin  corresponds  to  a current  source 
of  one -fourth  normal  value  if  only  one  quadrant  of  the  p- plane  is  represented 
by  the  analog.  IfV2tte  p-plane  is  used,  the  calibrating  singularity  will  have 
half  the  normal  value  of  current. 
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No  mention  has,  as  yet,  been  made  of  setting  the  vertical  control. 
It  will  be  recalled  that  one  of  the  variables  of  the  approximations  is  the 
determination  of  the  multiplicative  constant  of  the  imraittance  (Eq.  (l)). 
When  the  logarithm  of  the  amplitude  is  taken,  this  becomes  an  additive  con- 
stant (Eq.  (2)).  But  the  addition  to  a constant  of  a curve  means  simply 
raising  or  lowering  the  position  of  the  curve.  Thus,  the  vertical  position 
control  is  a multiplicative  adjustment  in  the  approximation  procedure. 


IX-  The  Approximation  from  the  Viewpoint  of  £conoml2ati on 

If  the  approximation  has  its  start  in  an  arbitrary  choice  of  the 
number  of  poles  and  zeros  and  a random  location  of  these  critical  points,  the 
problem  becomes  a question  of  relocation  of  these  points  to  obtain  the  best 
approximation  and  the  convergence  technique  is  not  readily  apparent.  In  fact, 
the  solution  of  this  problem  involves  so  much  ingenuity  that  it  may  at  times 
seem  to  be  more  of  an  art  than  a science.  For  this  reason,  the  approxima- 
tions carried  out  in  this  paper  start  with  a known  approximating  function.  An 
economization  process  is  used  to  reduce  the  number  of  poles  and  zeros  as  much 
as  possible,  thus  minimizing  the  number  of  elements  and  obtaining  the  best  ap- 
proximation for  this  situation.  The  method  uses  several  of  the  critical  points 
of  the  known  approximating  function  as  a starting  point.  Enough  poles  and 
zeros  are  chosen  to  give  a good  approximation  in  the  range  of  interest.  Then, 
the  pole  or  zero  having  the  least  effect  on  the  response  is  removed  and  the  re- 
maining critical  points  are  perturbed  until  the  response  is  returned,  as 
closely  as  possible.,  to  that  which  existed  before  the  removal  of  the  singularity. 
The  process  of  removal  of  one,  or  a pair,  of  critical  points  and  perturbation 
of  the  remaining  ones  is  then  repeated.  Usually,  the  amount  of  rearrangement 
necessary  after  the  removal  of  a pole  or  zero  is  not  too  great  until  a minimum 
number  of  critical  points  i3  reached.  Further  economization  makes  approximation 
within  specifications  impossible.  At  this  point,  the  optimum  synthesis  is 
reached.  Analytic  economization,  sometimes  called  "quantization",  has  also  been 
used  (References  5 and  12)  to  approximate  continuous  pole  and  zero  distributions 
with  a finite  number  of  critical  points. 


X.  The  Perturbation  Technique 

Although  it  is  not  always  possible  to  predict  the  movement  of  criti- 
cal points  in  order  to  converge  to  the  required  approximation,  a limited  per- 
turbation technique  has  been  evolved.  An  analytic  approximation  technique, 
by  successive  adjustments  specified  in  terms  of  Clifts  of  pole  and  zero  posi- 
tions made  to  reduce  deviation  from  desired  characteristics,  has  been  success- 
fully carried  out  elsewhere  (Reference  7).  However,  the  methods  are  much  more 
laborious  than  those  described  here  and  a check  on  optimum  synthesis  cannot 
readily  be  made. 
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The  first  question  in  the  economization  concerns  itself  with  the 
order  in  which  the  poles  and  zeros  are  removed.  In  general,  the  critical 
points  are  eliminated  as  follows: 


1.  Points  outside  the  region  of  interest. 

2.  Points  lying  at  the  far  ends  of  the  region  of  interest. 


3-  Points  furthest  from  the  real  frequency  axis. 

U.  Current  sources  of  one  polarity  lying  between  sources  of  opposite 

polarity. 

5.  Other  points,  if  possible. 

Step  U is  perhaps  most  in  need  of  explanation.  It  sometimes  occurs 
that  three  critical  points,  say  a zero  between  two  poles,  lie  along  a line 
approximately  parallel  to  the  scanning  axis  (constant  phase  line).  It  is  then 
possible  to  economize  by  removing  the  center  critical  point.  The  success  of 
the  economization  depends  somewhat  on  the  original  positions  of  the  critical 
points  and  on  the  total  numbers  of  poles  and  zeros  in  the  finite  region  of  the 
plane.  Another  point  of  interest  is  the  method  of  perturbation  which  leads  to 
a most  rapid  convergence.  An  empirical  technique,  applicable  in  many  situa- 
tions, has  been  evolved.  It  should  be  said  here  that  control  of  the  multiplica- 
tive constant  plays  an  important  role  in  the  convergence  procedure.  Vertical 
positioning  is  adjusted  following  every  other  step,  that  is,  after  each  removal 
of  a pole  or  zero;  again  after  each  perturbation  of  each  remaining  critical 
point. 

The  critical  points  may  be  classified  according  to  their  location  in 
the  logarithmic  plane  as: 


1.  Poles  or  zeros  close  to  the  real  frequency  axis.  These  points  lie 
close  to  maxima  or  minima  of  the  required  curve  or  where  the  curve  has  a steep 
slope.  The  position  of  these  points  is  most  critical,  therefore  they  are  the 
first  to  be  adjusted.  Moving  a pole  closer  to  the  scanning  axis  (along  a con- 
stant amplitude  line)  increases  the  peak.  Moving  the  pole  parallel  to  the  scan- 
ning axis  (along  a constant  phase  line)  keeps  the  amplitude  of  the  peak  essen- 
tially constant  but  changes  the  frequency  at  which  the  peak  occurs.  The  effect 
of  small  movement  along  a constant  amplitude  line  is  not  as  readily  observable 
far  from  the  critical  frequency  as  is  small  movement  along  a constant  phase 
line. 


2.  Poles  or  zeros  far  from  the  real  frequency  axis.  The  points  close 
to  the  axis  located  n/2  radians  from  the  scanning  line,  have  little  effect 
on  the  peaks  and  troughs  of  the  response  curve.  They  do  have  an  effect  over 
the  region  of  frequencies  nearest  to  their  location.  Their  movement  generally 
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raises  or  lowers  that  part  of  the  curve  where  the  slope  is  not  very  steep  and 
does  not  change  rapidly. 


3.  The  effect  of  movement  of  critical  points  centrally  located  in  the 
logarithmic  plane  is  most  difficult  to  predict  and  is  most  dependent  on  the 
particular  distribution  in  question. 

Two  general  statements  can  be  now  made.  Unless  there  is  great  con- 
centration or  bunching  of  poles  and  zeros,  then  for  each  critical  point  there 
is  a finite  range  of  frequencies  where  the  curve  is  critically  affected  by  a 
small  movement  of  the  point.  The  movement  is  small  enough  so  that  it  does  not 
greatly  affect  the  response  outside  this  finite  range.  Furthermore,  the  value 
of  the  response  can  be  kept  constant  at  any  frequency  even  if  all  poles  and 
zeros  are  moved,  provided  that  the  distances  to  these  critical  points  from  the 
frequency  in  question  is  kept  constant. 

Keeping  these  facts  in  mind,  it  becomes  a matter  of  experience  to 
obtain  the  skill  needed  for  making  the  perturbations  so  as  to  obtain  rapid 
convergence . 


n.  Synthesis  of  Driving-Point  Functions 

A necessary  and  sufficient  condition  for  the  physical  realizability 
of  a driving-point  Irani ttance  is  that  it  be  a positive  real  function.  This 
condition  is  guaranteed  if  the  approximation  is  carried  out  in  the  real  part 
plane  (to  b9  described  below),  providing  certain  restrictions  on  the  positions 
of  the  critical  points  are  satisfied. 


The  input  immittance  may  be  written  as 

»o  ♦ *iP  + •••  + aipl 
bQ  + b^  + . . . + bjp^ 

mg  + Joxig 

where  ra  and  n are  polynomials  even  in  co.  Then, 

2 

ro^nig  co 

— 2 — T2 — 

nig  + co  ng 


F(p)  - 

Along  jco  this  becomes 

F(co)  - 


(17) 


(18) 


Re  F(co) 


(19) 
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is  an  even  function  of  <0.  Substituting  X for  « gives 
Re  F(\)  - 


2 2i 

Aq  + AgX  ♦ . . . ♦ Ag^X 

B0  ♦ BgX2  ♦ ...  ♦ B2JX2j 


Consider  X to  be  a complex  variable 


X - co  ♦ J-< 


and  express  Re  F(X)  in  terms  of  its  critical  points 


Re  F(  X)  - H 


*(X2  - X2) 

*(x2  - A*) 


where  XQ  and  X^  are  the  zeros  and  poles,  respectively.  Factoring 

n(X  - \)(\  * XQ) 


Re  F(X)  - H 


"(a  - ycx  ♦ Xp) 


(20) 


(21) 


(22) 


(23) 


But  the  coefficients  of  Eq.  (20)  are  real  so  that  the  poles  and  zeros  of  the 
function  are  complex  conjugate  or  real.  Then  terms  containing  the  conjugate 
critical  points  of  Eq.  (23)  must  appear  whenever  these  points  are  complex,  or 


Re  F(X) 


(A  - XQ)(X  - A*)(A  + XQ)(X+  X*) 
(A  - Ap)( A - x£)(A  ♦ Xp) (A  ♦ Xj) 


(2U) 


Thus,  it  is  seen  that  the  poles  and  zeros  are  symmetrically  positioned  with 
respect  to  the  real  and  imaginary  axis  in  the  X- plane.  In  order  that  the  real 
part  be  associated  with  a physical  iranittance  function,  there  are  two  re- 
strictions made  on  the  locations  of  the  poles  and  zeros. 

1.  There  may  be  no  poles  anywhere  on  the  real  X axis  over  the  whole 
frequency  spectrum. 

2.  Zeros  on  the  real  axis  must  be  of  even  multiplicity.  Critical  points 
may  appear  anywhere  else  provided  that  the  symmetry  conditions  are  satisfied. 

Symmetry  with  respect  to  both  the  real  and  the  imaginary  axis  enables 
the  analog  to  represent  the  whole  function  with  any  quadrant  of  the  X-plane,  as 
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explained  in  Chapter  IV.  To  assure  physical  realizability,  the  approximation 
is  made  to  the  real  part  of  the  required  function.  The  necessary  restrictions 
on  the  locations  of  the  poles  and  zeros  are  very  easily  met  in  the  analog.  No 
pole-currents  are  introduced  on  the  scanning  axis  and  zero-currents  on  co  are 
even  multiples  of  unit  current.  In  order  that  no  pole  appear  at  infinity,  no 
pole-currents  may  be  introduced  at  the  infinity  return  equipotential.  This 
means  that  the  number  of  zeros  in  the  infinite  region  must  be  at  least  equal  to 
the  number  of  poles. 

Simple  analytic  methods  for  determining  F(p)  as  a minimum  reactance 
function  are  known  (Reference  9)  once  the  equation  of  Re  F(X)  has  been  found. 


XII.  The  Experimental  Results 

The  input  admittance  of  a short  circuited  transmission  line  was  syn- 
thesized over  a specified  frequency  range  with  a lumped-parameter  network  by 
approximating  to  the  input  conductance. 

The  input  admittance  of  a short  circuited  transmission  line  is  given 
by 


Y(oj)  - /..G  — ^°3G.  Coth  JL  J(R  - JML)  (G  - joC)  (25) 

V R + joiL 


where  R,  G,  C,  and  L are  the  resistance,  conductance,  capacitance,  and  induc- 
tance per  unit  length.  By  neglecting  the  conductance,  introducing 


Z 


c 


J L/C  lossless  characteristic  impedance 


f 

r 


- RjP/Z  attenuation  parameter 

c 

- 2n$/X  « cdJZ/v  electrical  argument 

- 1/  y/  LC  velocity  of  propagation 


normalizing  with  respect  to  Z and  substituting  s for  j T , Eq.  (25)  is  con- 
verted into 

Coth  s / 1 + -^/s 
/ (1  ♦f/s) 


Ks) 


(26) 
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The  Mittag-Leffler's  partial  fractions  expansion  of  Y(s)  is  (Reference  11) 


1(8) 


2s 

”2 TT 

s ♦ ^ s ♦ n n 


(27) 


A.  Approximation  for  the  Range 3.0 


The  approximation  in  this  range  was  achieved  by  means  of  an  econo- 
mization process  starting  with  the  first  two  terms  of  the  partial  fractions 
expansion.  The  starting  approximating  function  is 

I (s)  - — + 2s (28) 

^ ♦ 8 8 +^S  + It 

for  real  frequencies  this  becomes. 


Y (T)  - 1-08$  (9.87  - 3r2  ♦ 1 1-3U3 

1 1*.1*1  -.891*1*  ♦ j 1(10.07  -T2) 


(29) 


where  the  multiplicative  constant  1.085  was  chosen  in  the  analog  so  as  to  make 
Ox(Tr)  an  equal  ripple  approximation  to  Re  Y('s)  for'fc'S.  3.0. 


s 

R.  hm  - ffl(r)  . 1.08S  C1-3U-^  -B.i 1 t J (30) 

(1*.1*1  - .891*1*)  *'Z2  (10.07  -fc2) 

Substituting  X for  T and  letting  X take  complex  values,  the  critical  points  of 
G1(X)  located  in  the  first  quadrant  arej 

zeros:  Xg  - 2.39  | 2$° 

X1  - .1*1*7  I 90° 

X3  - 3.1 1*  I 1*.08° 


poles: 


(31) 
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These,  then,  were  the  original  approximating  critical  points.  The 
resultant  c&n«kXc£  function  is  to  approximate  Re  T(s)  f or  T <L  3.0  in  an  equal 
ripple  manner.  The  curve  to  which  the  approximation  was  made  is  a plot  of 
the  real  part  of  the  function  T(s)  of  Eq.  (26).  In  this  synthesis,  the  econo- 
misation led  to  a very  convenient  form.  It  was  found  that  the  approximation 
to  the  real  part  could  be  made  with  a function  having  poles  in  the  finite  re- 
gion and  with  all  the  seres  moved  to  infinity.  This  configuration  of  singulari- 
ties is  advantageous  since  the  resulting  network  can  always  be  synthesized  with 
a reactive  ladder  structure  containing  no  coupled  coils  and  terminated  in  a 
pure  resistance. 

Consider  Fig.  1. 


FIGURE  1 


The  short  circuit  transfer  admittance  iB 

712 


12 


then. 


'12 


y22  + 1 
N 

®2  * Pn2 

N/pn2 

mg/pag  ♦ 1 

N 

Fn2 


(32) 


(33) 


R-327-53,  PIB-263 


18 


and. 


'22 


is  purely  reactive. 


(3U) 


pnc 


Now,  the  zeros  of  power  transfer  are  also  the  zeros  of  sins 
at  the  input  terminals  1,  results  in  an  input  power  of 


v ClvdgO  f 


Ivl2  0, 


- hnr,jz 


'in  " '*12 

when  terminals  2 are  shorted.  Or, 

„ i2 


(35) 


in 


l12 1 


(36) 


If  all  the  zeros  of  Gin  are  at  infinity,  then  N is  a constant. 


Darlington  synthesis  (Reference  13)  leads  to  an  expansion  of  yz2  in 
a continued  fraction  such  that  the  zeros  of  yi2  are  all  at  infinity.  The  re- 
sult is  the  required  network,  within  a multiplicative  constant,  in  the  form 
of  a reactive  four- pole  terminated  in  a resistor. 


The  results  of  the  approximation  are  shown  in  Figs.  MKI-13299  ard 
MKE-13300  for  C-  and  I,  respectively.  It  is  noted  that  two  terms  of  the  partial 
fractions  expansion  (with  the  proper  multiplicative  constant)  is  a good  ap- 
proximation to  the  required  response.  Removal  of  the  zero,  without  perturba- 
tion of  the  poles  results  in  a very  large  deviation  of  Gi  in  the  region 
2.5><2<3»0*  The  final  approximation  with  the  analog  computer  generally  gives 
better  results  than  the  partial  fractions  expansion  in  the  given  range  and  the 
resulting  network  contains  fewer  elements.  The  networks  arising  from  the  two 
approximations  are  given  in  Fig.  MRI-13301*  the  equation  of  the  aialog  approxi- 
mation being 


Y2(s) 


/ 

p 

s ♦ ,62s  ♦ 3.1U 

.36Us3  ♦ ,226s^  ♦ 3.U8s  ♦ 1.U5 


(37) 


A comparison  of  the  critical  points  in  the  two  approximations  is  given  in 
Table  1.  The  response  curves  of  log  amplitude  against  log  frequency,  as  seen 
on  the  cathode  ray  tube,  are  shown  in  Fig.  MRI-13297. 
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B-  Approximation  for  the  Range  '£'£6.0 

Having  arrived  at  satisfactory  results  for  the  range  3. 0,  the 
next  step  was  an  extended  range  $ the  range  of  approximation  was  doubled.  The 
addition  of  another  term  of  the  partial  fractions  expansion  gave  the  starting 
point  approximation  function, 


*3(s) 


_J 

^ ♦ s 


♦ 


2s 

. 2 

a a * n 


♦ 


2s 

|/s  + 1*jT 


(38) 


Choosing  a multiplicative  constant  of  1.31  so  that  the  poles  of  the  three 
terms  of  Y3(s)  give  an  equal  ripple  approximation  to  Re  Y(s)  forTf£;6.0, 
Eq.  (38)  becomes. 


Y3(s) 


1.31  (5s1*  * U.U7s3  ♦ lU9.1s2  * 66.18  ♦ 390) 
s5  ♦ 1.3UlsU  ♦ U9.97s^  ♦ UU.2s^  ♦ 399  s ♦ 17U 


f 


(39) 

.UU7 


The  desired  analog  function  is  to  approximate  Re  Y(s)  for  6.0 

in  an  equal  ripple  manner  with  all  zeros  in  the  real  part  plane  located  at 
infinity.  The  poles  of  Re  Y3(s)  give  the  unperturbed  location  of  the  singu- 
larities. These  are  obtainable  from  the  solutions  of  the  equations 


From  which  the  poles 


^ * jX  - 0 


n2  - X2  ♦ 

JfX  - 0 

m 

Un2  - X2  ♦ 

1 

0 

in  the  first 

quadrant  of  X are: 

\ - .1*2*7 

| 90° 

x3  - 3.1U 

| U.08° 

(U 1) 

- 6.28 

1 2.0U° 

The  results  of  the  analog  approximation  are  shown  in  Figs.  MRI-13302 
and  MRI °13303  for  0 and  Y,  respectively.  Fig.  MRI-13302  clearly  shows  that 
moving  the  zeros  of  the  partial  fractions  approximation  to  infinity  results  in 
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large  deviations  from  the  required  response.  Adjustments  of  pole  positions, 
however,  result  in  an  analog  approximation  which  is  considerably  better.  Fig* 
MRI-I33O3  gives  the  resultant  plots  of  the  amplitude  of  the  admittance.  Al- 
though the  analog  approximation  deviates  from  the  exact  amplitude  of  the  ad- 
mittance to  a marked  degree  in  the  vicinity  of  the  peak  (at't  ■ 3.H*)»  it  is 
considerably  better  than  the  partial  fractions  expansion.  Furthermore,  an 
economization  of  two  elements  has  been  affected,  and  the  analog  network  is  of 
simpler  configuration.  The  networks  arising  from  the  two  approximations  are 
given  in  Fig.  MRI-1330ii,  the  equation  of  the  analog  approximation  being 


V8> " 


.191s 


j1*  ♦ .8?8s3  + 36s2  ♦ 2U.7s  ♦ 76.U 


1673^  ♦ 8.9iia3  + 6.5Us2  ♦ 68. Us  ♦ 33^2 


(U2) 


It  should  be  emphasized  that  the  analog  approximation  is  made  to  the  real  part, 
and  the  resultant  approximation  to  the  amplitude  is  computed  by  constructing 
the  complete  complex  admittance  according  to  References  9 and  13. 

A comparison  of  the  critical  points  in  the  two  approximations  for 
the  extended  range  is  given  in  Table  2.  It  is  to  be  observed  that  relatively 
small  adjustments  were  required  to  obtain  the  improvements.  The  response 
curves  of  log  amplitude  against  log  frequency,  as  seen  on  the  cathode  ray 
tube,  are  shown  in  Fig.  MRI-13298. 


HII.  Conclusion 


An  electrostatic  analog  computer,  capable  of  approximating  to  the 
real  or  imaginary  part  as  well  as  to  the  amplitude,  has  been  successfully 
built.  The  machine  gives  results  well  within  experimental  accuracy  and  is  a 
valuable  laboratory  tool.  Although  the  method  is  one  of  trial- and- err or,  the 
amount  of  computational  work,  as  compared  to  that  normally  encountered  in  the 
analytical  approach  to  the  approximation  problem,  is  greatly  reduced. 

The  perturbation  technique,  as  well  as  the  proper  use  of  the  log- 
arithmic transformation,  simplifies  the  approach.  The  continuous  scan  feature 
of  the  analog  makes  it  possible  to  immediately  observe  changes  in  the  approxi- 
mation and  solutions  are  obtainable  in  a very  short  time.  The  approach  from 
the  real  part  leads  to  realizable  driving-point  functions,  while  economization 
reduces  the  number  of  elements.  The  experimental  results  show  that  relatively 
small  adjustments  in  the  positions  of  the  critical  points  affected  substantial 
improvements.  It  should  be  pointed  out  that  the  same  approximations  to  the 
transmission  line  could  have  been  gotten  without  a knowledge  of  the  partial 
fractions  approximating  networks  since  the  location  of  critical  points  are 
readily  obtainable. 
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